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C. Schwanda* (for the Belle collaboration) 

High Energy Accelerator Research Organization (KEK), Tsukuba, Japan 

We have searched for the decay B + — > ujI + v in 78 fb _1 of T(45) data (85.0 million BB events) accumulated with the 
Belle detector. The final state is fully reconstructed using the uj decay into ir + ir~ir° and the detector hermeticity to 
infer the neutrino momentum. The signal yield is extracted by a two-dimensional fit to the lepton momentum and the 
invariant 7r + 7r _ 7r° mass. The result of the fit depends on the form factor model assumed for the decay. Taking the average 
over three different models, 421 ± 132 events are found in the data, corresponding to a preliminary branching fraction of 
(1.4 ± OA(stat) ± 0.2(syst) ± 0.3(model)) ■ 10 -4 . 



1 Introduction 

Tree-level semileptonic B decays are crucial for mea- 
suring the CKM matrix elements \V c b\ and \V u b\ [H] 
independently of any new physics contribution. The 
decay B — > X u lv 1 (which allows to access \V u b\) has 
been studied both using inclusive approaches (sensi- 
tive to all X u lv final states within a given region of 
phase space) and through exclusive reconstruction of 
specific final states [E]. As to the latter, Belle has 
already obtained preliminary results for the decays 
B° -> Ti-l+v and B+ -> p°l+v [ 3 . In this article, we 
present evidence for the decay B + — ► ujI + v which is 
so far unobserved. 



2 Experimental procedure 

2.1 KEKB and the Belle detector 

Belle is located at the KEKB asymmetric e + e~ col- 
lider, operating at the T(AS) resonance [ 0J. The 
detector is described in detail elsewhere [ It is 
a large-solid-angle magnetic spectrometer that con- 
sists of a three-layer silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T mag- 
netic field. 

The responses of the ECL, CDC (dE/dx) and ACC 
detectors are combined to provide clean electron iden- 
tification. Muons are identified in the instrumented 
iron flux-return (KLM) located outside of the coil. 



Charged hadron identification relies on information 
from the CDC, ACC and TOF. 

2.2 The data set 

This analysis is based on an integrated luminosity of 
78.13 fb _1 taken on the T(4S I ) resonance. This data 
set contains (85.0 ±0.5) • 10 6 BB events. Furthermore, 
to subtract the hadronic background, 8.83 fb _1 taken 
below the resonance are used. 

Full detector simulation is applied to Monte Carlo 
events. Generic background Monte Carlo samples, 
equivalent to about three times the beam luminos- 
ity, are used in various stages of this analysis. Monte 
Carlo samples for the signal, B + — > loI + v, are gener- 
ated with the following form factor models: ISGW2 
(quark model [©), UKQCD (quenched lattice QCD 
calculation [[7]) and LCSR (light cone sum rules 
The cross-feed from other decays B — ► X u lv is esti- 
mated using the ISGW2 Monte Carlo. 

2.3 Neutrino reconstruction 

Events passing the hadronic selection are required to 
contain a single lepton (electron or muon) with a cm. 
momentum p\ greater than 1.3 GeV/c. In this mo- 
mentum range, electrons (muons) are selected with 
an efficiency of 92% (89%) and a pion fake rate of 
0.25% (1.4%). This requirement is 41% efficient for 
B — > X u lv events. 

The missing four-momentum is computed for selected 
events, 



= Pher + Pler 



* Supported by the Japan Society for the Promotion of Science. 
1 Throughout this article, the inclusion of the charge conjugate 
mode is implied. 



E m iss — Eher + Eler — ^ Ei , (1) 



where the sum runs over all reconstructed charged 
tracks and photons. The indices HER and LER re- 
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Figure 1. The missing momentum resolution in B — > 
X u li/ events (ISGW2 model). The fit shown uses a double 
Gaussian (<xi = 141 ± 10 MeV/c, <t 2 = 412 ± 8 MeV/c). 



fer to the high energy and the low energy ring, re- 
spectively. To reject events in which the missing mo- 
mentum does not match the neutrino momentum, the 
following selection criteria are applied. Events with 
a large charge imbalance are eliminated, |Qtot| < 3e, 
and the direction of the missing momentum is required 
to lie within the ECL acceptance, 17° < 9 m iss < 150°. 



miss miss 



Pr, 



IS 



The missing mass squared, m 
required to be consistent with the neutrino hypothesis, 
\m 2 miss \ < 3 GeV 2 /c 4 . 

After applying these cuts, the resolution in p m i SS 
(E miss ) is found to be 141 MeV/c (587 MeV) for sim- 
ulated B — » X u lv events (Fig.Q. As the energy res- 
olution is worse than the momentum resolution, the 
neutrino four- momentum is taken to be (j> m iss,Pmiss)- 
The efficiency of the event selection and neutrino re- 
construction cuts is about 17% for B — » X u lv events. 



2.4 Final state reconstruction 

Pairs of 7's are combined to form 7r° candidates (E^ > 
30 MeV, 120 < 771(77) < 150 MeV/c 2 ). The de- 
cay uj — » 7r + 7r~7r° (branching ratio: (89.1 ± 0.7)% [ 
El) is reconstructed by trying all possible combinations 
of one ir° with two oppositely charged tracks. Com- 
binations with a charged track identified as a kaon 
are rejected and the following selections are imposed: 
p* > 300 MeV/c, 703 < m(7r+7r-7r°) < 863 MeV/c 2 . 
Combinations located far away from the center of the 
Dalitz plot are removed by requiring the Dalitz am- 
plitude, A oc \p^+ x fin- I, to be larger than half of its 



maximum value. 

The single lepton in the event is combined with the 
to candidate and the inclusive neutrino, and the lepton 
momentum cut is tightened, 1.5 < p ; * < 2.7 GeV/c. 
To reject combinations inconsistent with signal decay 
kinematics, the cut |cos#sy| < 1.1 is imposed, 



cos 9 



BY 



9 Z?* Z7 1 * 



2p* B p* Y 



(2) 



where Eg, p* B and uib are fixed to their nominal val- 
ues, VEherEler, v^b " m B and 5 - 279 GeV/c 2 , 
respectively. Ey, p Y and my are the measured cm. 
energy, momentum and mass of the Y = lu + 1 system, 
respectively. For well-reconstructed signal events, 
cos 9 by is the cosine of the angle between the B and 
the Y system and lies between —1 and +1 while for the 
various backgrounds, a significant fraction is outside 
this interval. 

For each ujIv candidate, the beam-constrained mass 
m bc and AE are calculated (££ eom = VEherEler), 



^(ELam) 2 -\Pl+Pl+Pl\ 2 



AE = Ei 



(e:+e; + e:) 



(3) 



and candidates in the window m{, c > 5.23 GeV/c 2 and 
\AE\ < 0.36 GeV are selected. On the average, two 
combinations per event remain after applying all cuts. 
We choose the one with the largest ui momentum in 
the cm. frame, which is the right choice about 80% of 
the time. 

2.5 Continuum suppression 

The hadronic continuum from e + e~ — > qq events, 
q = u,d,s,c, is suppressed using three variables ex- 
ploiting the fact that, in the T(45) frame, the two 
B mesons are produced nearly at rest and that there- 
fore BB events have a nearly spherical shape while 
continuum events have a more jet-like topology. These 
variables are: 

• The ratio R2 of the second to the zeroth Fox- 
Wolfram moment [ [23 . This ratio tends to 
be close to zero (unity) for spherical (jet-like) 
events. (The cut R2 < 0.4 is imposed already at 
the event selection level.) 

• The cosine of 9 thrust, where 9 t hrust is the angle 
between the thrust axis of the ljI system and the 
thrust axis of the rest of the event. 

• A Fisher discriminant selecting events with an 
even energy distribution around the lepton di- 
rection [^2- The input variables are the charged 
and neutral energy in nine cones of equal solid 
angle around the lepton momentum axis. 
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To optimize their individual discriminating power, 
these three variables are combined into a likelihood 
ratio. This selection is 55% efficient for B + — > ul + v 
events, while 92% of the continuum background re- 
maining after the R2 < 0.4 cut, is eliminated. 

2.6 The fit 

The amount of signal and the remaining background 
in the mf, c vs. AE signal window are determined by 
a two-dimensional binned likelihood fit to the lepton 
momentum p\ (4 bins, width: 300 MeV/c) and the in- 
variant mass m(7T + 7r~7r°) (8 bins, width: 20 MeV/c 2 ) [ 
112) . The following five components are fitted to the 
data: B + — > uil + v signal, B — > X u lv background, 
B — > X c lv background, fake and non B decay lepton 
background and continuum. The shapes of the for- 
mer four components are determined from the simula- 
tion, the shape of the continuum component is given 
by the off-resonance data. The normalizations of the 
B + — > ujI + v, the B — » X u lv and the B — ► X c lv com- 
ponent are floated in the fit, the other components 
being fixed (Fig. 0) . 

3 Result and systematic uncertainty 

The signal yield, N(B + — > ujI + v), is determined sep- 
arately for each form factor model and the branching 
ratio, B(B + — > iol + v), is calculated using the relation 

N{B + -> wl + v) = N(B+) x B{B+ -> ujI + v)x 

xB{u ^ ir+n-ir ) x (e e + e p ) , (4) 

where N(B + ) is the total number of charged B mesons 
in the data, B(u -> tt + tt~tt°) = (89.1 ±0.7)% [El and 
e e (e M ) is the model dependent selection efficiency for 
the electron (muon) channel (Tabled . Averaging over 
the three models (giving equal weight to each), 421 ± 
132 signal events are found in the data, corresponding 
to a branching fraction of (1.4 ± OA(stat)) ■ 10" 4 . The 
spread between the different models amounts to 21% 
and is used as an estimate of the form factor model 
uncertainty. 

The largest experimental uncertainty is the efficiency 
of the neutrino reconstruction (7.5%). It is estimated 
by varying the track finding efficiency, the neutral 
cluster finding efficiency, the track momentum res- 
olution and the cluster energy resolution separately 
in the Monte Carlo within their estimated range of 
uncertainty. The next to largest component is the 
cross-feed from decays B — > X u lv (4.4%) which is 
estimated by varying the fraction of B — > -nlv and 
B — > plv (these decays are expected to dominate in 
the high p[ region) within their experimental uncer- 
tainty. Other contributions are: B — > X c lv cross-feed 



(2.4%), charged track and cluster finding (3.0% and 
4.0%, respectively), lepton identification (3.0%), the 
number of BB events (0.6%) and the uncertainty in 
the u> — > 7T + 7r~7r° branching fraction (0.8%). 

4 Conclusion 

We have measured the decay B + — > ul + v using 
78 ib" 1 of T(45) data (85.0 million BB events). The 
final state of the decay was fully reconstructed using 
the oj decay into 7r + 7r-7r and detector hermeticity to 
infer the neutrino momentum. The signal yield was 
determined by a two-dimensional fit to the lepton mo- 
mentum and the invariant ir + ir~ir° mass. Repeating 
the fit for three different B + — > loI + v form factor mod- 
els and averaging the result, 421 ± 132 signal events 
are found, corresponding to a branching fraction of 
(lA±0A(stat)±0.2(syst)±0.3(model))-W- i for this 
decay. This result is preliminary. 
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Figure 2. The result of the fit in bins of p* , using ISGW2 form factors for the decay B + — » ajZ + ZA The data points are 
continuum subtracted on-resonance data, the histograms are the components of the fit, as described in the text. 



form factor model 


N(B+ -► ujI+v) 


B(B+ -> ul+v) 


x 2 


ISGW2 
UKQCD 
LCSR 


359 ± 106 
428 ± 134 
476 ± 156 


(1.02 ±0.30) • 10- 4 
(1.43 ±0.45) • 10" 4 
(1.73 ±0.57) • 10- 4 


30.5 
30.6 
30.6 


average 


421 ± 132 


(1.39 ±0.44) • 10~ 4 





Table 1. The signal yield, the corresponding B + — > ul + u branching fraction and the x 2 of the fit (ndf = 31), for each of 
the three form factor models used. 



